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Abstract 

Observation of top quark flavor changing neutral process t — > c + 7 at the LHC would be 
the signal of physics beyond the Standard Model. If anomalous t — > cj coupling exists, it will 
affect the precisely measured B{B — > X s ^f). In this paper, we study the effects of a dimension 
5 anomalous tcj operator in B — > Xg'j decay to derive constraints on its possible strength. It 
is found that, for real anomalous t — > 07 coupling k^ cR , the constraints correspond to the upper 
bounds B(t -> c + 7) < 6.54 x lO" 5 (for nJ cR > 0) and B(t -> c + 7) < 8.52 x lO" 5 (for ^ cR < 0), 
respectively, which are about the same order as the 5a discovery potential of ATLAS (9.4 x 10~ 5 ) 
and slightly lower than that of CMS (4.1 x 10~ 4 ) with 10 fb _1 integrated luminosity operating at 
y/s = 14 TeV. 
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I. INTRODUCTION 



In the Standard Model (SM), top quark lifetime is dominated by the t — > bW + process, 
and its flavor changing neutral current (FCNC) processes t — > qV(q = u,c;V = j,Z,g) 
are extremely suppressed by GIM mechanism. It is known that the SM predicts very tiny 
top FCNC branching ratio B(t — > qV), less than O(10~ 10 ) [lj], which would be inaccessible 
at the CERN Large Hadron Collider (LHC). In the literature as 

, however, a number of 

interesting questions have been intrigued by the large top quark mass which is close to the 
scale of electroweak symmetry breaking. For example, one may raise the question whether 
new physics (NP) beyond the SM could manifest itself in nonstandard couplings of top quark 
which would show up as anomalies in the top quark productions and decays. 

At present, the direct constraints on Bit — > qV) are still very weak. For its radiative 
decay, the available experimental bounds are Bit — > uj) < 0.75% from ZEUS [4J and 
Bit — > qj) < 3.2% from CDF at 95% C.L., respectively. These constraints will be 
improved greatly by the large top quark sample to be available at the LHC, which is expected 
to produce 8 x 10 6 top quark pairs and another few million single top quarks per year at 
low luminosity (10 fb~ 1 /year). Both ATLAS [6] and CMS {?| have got analyses ready for 
hunting out top quark FCNC processes as powerful probes for NP. With 10 fb _1 data, it is 
expected that both ATLAS and CMS could observe t — y qj decays if their branching ratios 
are enhanced to C(10~ 4 ) by anomalous top quark couplings (], 7|. However, if the top quark 
anomalous couplings present, they will affect some precisely measured qualities with virtual 
top quark contribution. Inversely, these qualities can also restrict the possible number of top 
quark FCNC decay signals at the LHC. The precisely measured inclusive decay B — > X s ^ 
is one of the well known sensitive probes for extensions of the SM, especially the NPs which 
alter the strength of FCNCs s| . Thus, when performing the study of the possible strength 
of t — > cy decays at the LHC, one should take into account the constraints from B — > 

In this paper, we will study the contribution of anomalous t^yc operators to the B — > X s 7 
branching ratio and derive constraints on its strength. In the next section, after a brief 
discussion of a set of model-independent dimension 5 effective operators relevant to t — > 07 
decay, we calculate the effects of operator c^o^t^F^ in B — > X s 7 decay, which result in a 
modification to CV 7 . In Sec. Ill we present our numerical results of the constraints on its 
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strength and the corresponding upper limits on branching ratio of t — > cy decays. Finally, 
conclusions are made in Sec. IV. Calculation details are presented in Appendix A, and input 
parameters are collected in Appendix B. 



II. TOP QUARK ANOMALOUS COUPLINGS AND THEIR EFFECTS IN B 
X,7 DECAY 

Without resorting to the detailed flavor structure of a specific NP model, the Lagrangian 
describing the top quark anomalous couplings can be written in a model independent way 



with dimension 5 operators 



III 



K 9 n k w ft 7 



A ™ u v J2 ^ A w " ^ ^ A 

q=u,c,t v q=d,s,b q=u,c,t 



- E ^QR^lZ^ + (R^L) + h.c, (1) 

2 cos Vw A 

q=u,c,t 

where k is the complex coupling of its corresponding operator, 6\y is the weak angle, and 
T a is the Gell-Mann matrix. A is the possible new physics scale, which is unknown but may 
be much larger than the electroweak scale. There are also Lagrangian describing the top 
quark anomalous interactions with dimension 4 and 6 operators, and the dimension 4 and 
5 terms can be traced back to dimension 6 operators [3, Q|- In fact top quark anomalous 
interactions can be generally described by the gauge-invariant effective Lagrangian with 
dimension 6 operators in a form without redundant operators and parameters Q, 
A recent full list of dimension 6 operators could be found in Ref. |L5|] . But for on-shell 
gauge bosons, the Lagrangian in Eq. ([!]) works and is commonly employed in high energy 



phenomenology analysis [3 



The operators in Eq. ([I]) relevant to t — > q^y decays read 



7 7 
K+„ r x — -\ ft+ 



= ~ e E -f-^t^LF^ - e E -f-QL^tnF^ + h.c. (2) 

q=u,c q=u,c 

It is understood that the Dirac matrix a^ u connects left-handed fields to right-handed 
fields, the t cry transition will involve two independent operators Vfiqq^a^t^F^ and 
m t(lL( J ' J ' u tRF f j iU , where the mass factors must appear whenever a chirality flip L — > R or 
R — > L occurs. Due to the mass hierarchy m t ^> m c , the effect of m q q R a^ u t L F^ v can be 
neglected unless n] qL is enhanced to be comparable to ^k^r by unknown mechanism. 






q = u,c q = u, c 



(a) (b) (c) 

FIG. 1: Feynman diagrams for b — > sj. (a) and (b) are the penguin diagrams with the anomalous 
tqj couplings, (c) Sample LO penguin diagram in the SM. 

The anomalous t*yq coupling affects b — » S7 decays through the two Feynman diagrams 



depicted in Figs. I (a) and 1(b) It is interesting to note that the CKM factors in Fig. l(a 



and Fig. I (b) are V tb V* and V qb V* s , respectively. Since |V i6 V*| > (V^V^I for q = u,c, the 



contribution of Fig. 1(a) would be much stronger than that of Fig. 1(b) Furthermore, given 



the strengths of t — >■ 1*7 and t — > 07 comparable, the contribution of Fig. 1(a) to b — > 57 is 



still dominated by t — > cy because of |V CS | 3> \V US \. Hence we will only consider Fig. 1(a) with 



anomalous tcy coupling. From the Feynman diagram of Fig. 1(a), it is easy to observe that 
the large CKM factor V t bV cs ~ 1 makes b — > 57 very sensitive to the strength of anomalous 
tcy coupling. 



The calculation of Fig. 1(a) can be carried out straightforwardly. The calculation details 



are presented in Appendix A, and the final result reads 

iM(b^ sj) = s[eT v (k)}be v (k), 
G f 



eY u {p } k) = ie- 



V c * s V tb [ia^k^mJML + m b f R (x)R)] 



Usually m s term can be neglected, and the function /r(x) is calculated to be 

1 



(3) 



K 
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(x c - l)(x t - 1) (x c - l) 2 (x c - x t ) 



In x c + 



X7 



(x t - l) 2 (x c - x t ) 



In x t 



4) 

with x q = m^/mly. Now we are ready to incorporate the NP contribution into its SM 
counterpart for B — > X s 7 decay. 

In the SM, it is known that B — > A s 7 decay is governed by the effective Hamiltonian at 
scale /i = 0{m b ) [17] 

" 6 

J^Cii^Qii^ + Cr^Oj^ + Csgi^Osgifi) , (5) 



n eS (b s 7 ) = -^v*v tb 



V2 



i=i 
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where Cj(/x) are the Wilsion coefficients, Oj = i_6 are the effective four quark operators and 



O 



n 



;m b (s L a^b R )F^, O 
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For calculating B{B X s j), instead of the original Wision coefficients Cj, it is convenient 
to use the so called "effective coefficients" [18j 



; m b (s L a^T a b R )G« . 



(6) 



< )cff K) = V*CW"(M W ) + „S)c£> s M (M W ) + C^ M (M W ) £ M"S (7) 



i=i 



where r\ = a s (iiw)/ a s(^b) and 
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a, = ( 
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-0.6494 - 0.0380 - 0.0185 -0.0057 ), (8) 
M 5s ~M °- 4086 -°- 4230 - 0.8994 0.1456 ). (9) 



2l|. 



To the leading order approximation, the B(B — > X s ^) is proportional to \C^L (7n&)| 2 

In terms of the operator basis in Eq. (j^J), the contribution of the anomalous t — > C7 
couplings in Eq. (j3J) would result in the deviation of 



C 7J (M W ) C 77 (M W ) = C 77 (M W ) + C 7 N 7 P (M W 



(10) 



and C7I (Mp^) can be read from Eq. as 



^7 N 7 P ( M ^) 



A V£ 



m t 



+ 



(x c - l)(x t - 1) (a; c - l) 2 (x c - x t ) 



\ogx c - 
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In x t 



(Xi - l) 2 (x c - x t ) 

(11) 

From this equation, one can see that the NP contribution is suppressed by a factor of m t /A 
but enhanced by V cs /Vt s . 

Since NP contribution does not bring about any new operator, the renormalization group 
evolution of from M w to mj scale is just the same as the SM one in Eq. ([7j). For m t = 172 
GeV, m b = 4.67 GeV, a s (M z ) = 0.118 and A = 1 TeV, we have 



16 
7723 



C\f M {M w ) + C^{M w ) 



O . 14 



^)< )SM (Mv,) + Cf SM (M w )^/ ij 



i=l 



0.665 



C\f M {M w ) + C^ iNt -(M w ) + 0.093 C^ ,om {M w ) - 0.158 C^ u;D1Vi (M w ) 



,(0)NP, 



(0)SM, 



,(0)SM, 



0.665 [-0.189 + «2dR(— 1.092)] + 0.093 (-0.095) - 0.158. 



( 



In principle, C'f^(m b ) will receive corrections from anomalous t — )■ eg couplings in Eq. ([I]) 
which will cause a deviation to C$ SM (M W ). However, as shown by Eq. f lT2"j) , the coefficient 



r/23 of C)^(M W ) is about one order larger than §(7723 — 7723) of Cg g (My/). Given the 
relative strength of C ( ^ V (M w ) to C^ )SM (M W ) at 10% level, C' 7 f(m b ) will be shifted by only 
few percentage. For simplifying the numerical analysis, we would neglect the contribution 
of the anomalous t — > eg couplings. We also find that the operator qua^t^F^ contributes 
to B — > Xs7 only through the term m s ser MI ,(l — 75)6 as shown by Eq. (|SJ) and Eq. (1A7I) . 
Combined with the previous remarks on this operator, the effects of qRO^t^F^ could be 
safely neglected. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



The current average of experimental results of B(B — > X^) by Heavy Flavor Average 
Group is \]J^ 

B cxp (B X sl ) = (3.55 ± 0.24 ± 0.09) x 10" 4 . (13) 



On the theoretical side, the NLO calculation has been completed [20j, |2l|, and gives 



B(B X s7 ) = (3.57 ± 0.30) x 10" 



(14) 



The recent estimation at NNLO 



22| gives B(B X s7 ) = (3.15±0.23) x KT 4 , which is about 



la lower than the experimental average in Eq. (113]) . Thus the experimental measurement 
of B(B — > Xs'y) is in good agreement with the SM predictions with roughly 10% errors on 
each side. The agreement would provide strong constraints on the top quark anomalous 
interactions beyond the SM 



The decay amplitude of t — > cry has been calculated up to NLO 



161 ] . For a consistent 



treatment of the constraints from t — > cy and b — > sj decays, we use the NLO formulas 



in Ref. 



201 ] to calculate B(B — > Xg^y). The experimental inputs and main formulas are 



collected in Appendix B. For numerical analysis, we will use the notation k^ cR = l^pJe^R 



and set A = 1 TeV. 

At first, we analyze the dependence of J g SM + NP (^ _ > X^) on the new physics parameters 
|fi^ cR /A| and ^ cR , which is shown in Fig. |2j From the figure, one can find that the contribu- 
tion of anomalous t — > cy coupling is constructive to the SM one for 9^ cR G [—50°, 50°], thus 
B(B — > X s ry) is very sensitive to |^ cR |. However, when |6^ cR | e [80°, 130°], the sensitivity 
of B(B — > X s ry) to |fi^ cR | becomes weak. For |6^ cR | ~ 180°, the contribution of anomalous 
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FIG. 2: The contour-plot describes the dependence of B(B — > X s 7)(xl0 -4 ) on |k^ cR /A| and 6^ cR - 
The dashed lines correspond to the experimental center value of B(B — > X s j). 

t — t- C7 coupling is destructive to the SM one and there are two separated possible strengths 
for Kr/A|- 
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FIG. 3: The 95% C.L. upper bounds on anomalous coupling |«^ cR /A| as a function of #^ cR - The 
shadowed region is allowed by B exp (B — > X s ^) and the dash-line is the CDF [5J upper limit. 

The allowed region for the parameters |k^ cR /A| and 6^ cR under the constraints from 
B(B — > X s ^) at 95% C.L. is shown in Fig. [3J The corresponding 95% C.L. upper bound on 
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FIG. 4: B(t — > C7) as a function of #^ cR . The shadowed region is allowed by the combined 
constraints of B(B — > X s j) and CDF searching at 95% C.L. 

Bit — > cry) is shown in Fig. HJ 

Now we turn to discuss the our numerical results. From Eq. ffP2"]) . the explicit relation 
between the SM and the t — > cy coupling contributions is 

C' 7 f(m b ) = -0.293 - 0.726 < cR . (15) 

Obviously, when Re k^ cR > 0, the interference between them is constructive, and it turns to 
be destructive when 6^ cR > 90°. Thus the features of these constraints shown in Figs. [3] and 
H]for different 6^ cR are 

(i) the bound on |k^ cR /A| is very strong for 6^ cR G [—50°, 50°]. For 6^ cR ~ 0°, as shown 
in Fig. [3J we obtain the most restrictive upper bound |/^ cR /A| < 4.9 x 10~ 5 GeV -1 , 
which implies Bit — >■ cy) < 6.54 x 10~ 5 ; 

(ii) the bound on |/t^ cR /A| is rather weak for 6^ cR around 110°. For such a case, Re fi^ cR ) is 
destructive to the SM contribution as shown by Eq. f)15p . so, the allowed strength for 
the anomalous coupling is much larger than the one for real k^ cR . When |6^ cR | ~ 135° 
and |fi^ cR | ~ 0.571, Cj^(m b ) is almost imaginary since Re C^{m b ) « 0. Then the 
restriction on |k^ cR /A| is provided by the CDF search for Bit — > cy) 

(iii) as shown in Fig. [3J when 6^ cR ~ ±180°, there are two solutions for |^ cR /A|. The 
larger one |/^ cR /A| ~ 1.4 x 10 -3 GeV _1 (S2 column in Tabled]) corresponds to the 
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TABLE I: The 95% C.L. constraints on the anomalous t cry coupling by B(B — > X s y) and 
B(t — > cry) for some specific #^ cR values. 

<%R = 0° 6^ = ±180° SI g^ = ±180°S2 ^ = ±110° 

B(B^X„i) \*l R \< 0.049 |< cR |< 0.056 1.35 < \kJ cR \ < 1.45 |«£ R | < 0.55 

B(t -> cry) CDF bounds[5] |^ cR | < 1.09 \kJ cR \ < 1.09 |^ cR | < 1.09 |re£ R | < 1.09 

Combined bounds |k^ cR | < 0.049 |k^ cR | < 0.056 - |k^ cR | < 0.55 

B(t -> cry) < 6.54 x 10~ 5 < 8.52 x 10~ 5 - < 8.17 x 10~ 3 

situation that the sign of is flipped. However, it has been excluded by the CDF 
upper bound of Bit — » cy) < 0.032 [5]. The another solution (SI column in Table H]) 
|/«t cR /A| < 5.6 x 10~ 5 GeV -1 will result in the upper limit Bit ->■ cy) < 8.52 x 10~ 5 . 

Taking 9^ cR = 0°, ± 180° and ±110° as benchmarks, we summarize our numerical con- 
straints on kJ cR and their corresponding upper limits on Bit — > 07) in Table [H From the 
table, we can find that our indirect bound on real «^ cR is much stronger than the CDF 
direct bound. The corresponding upper limits on Bit — > 07) are about the same order as the 
ATLAS sensitivity B(t 07) > 9.4 x 10" 5 Q and CMS sensitivity B(t -»■ 07) > 4.1 x 10" 4 
[3] with an integrated luminosity of 10 ftT 1 of the LHC operating at y/s = 14 TeV joj. 



IV. CONCLUSIONS 

In this paper, starting with model independent dimension five anomalous tcy operators, 
we have studied their contributions to B{B — > Xgj). It is noted that the t — >■ C7 transition 
will involve two independent operators n^CLcr^tnF^ and K^CRa^tLF^. The first oper- 
ator will produce a left-handed photon in t — > cy decay, while the second one will produce 
a right-handed photon. It is found that B — > X s y is sensitive to the first operator, but not 
to the second one. 

For real kJ cR , the constraint on the presence of «^ cR C£<7 Miy t rF^ is very strong, which 
corresponds to the indirect upper limits Bit — > cy) < 6.54 x 10 -5 (for positive k^cr) and 
Bit — > ay) < 8.52 x 10~ 5 (for negative nJ cR ), respectively. These upper limits for Bit — > cy) 



IT] 



are close to the 5a discovery sensitivities of ATLAS jo] and slightly lower than that of CMS 



with 10 fb 1 integrated luminosity operating at y/s = 14 TeV. For nearly imaginary k 



1 

tcR' 



the constraints are rather weak since CV 7 in the SM is a real number. If Bit — > cy) were 
found to be of the order of O(l0~ 3 ) at the LHC in the future, it would imply the weak phase 
of k^ cR to be around ±100°. However, such a coupling might be ruled out by the other 
observable in B meson decays jsoT ]. 

In summary, we have studied the interesting interplay between the precise measurement 
of b — > S7 decay at B factories and the possible t —> cry decay at the LHC. For real anomalous 
coupling, it is shown that B(t — > cy) has been restricted to be blow 10~ 4 at 95% C.L. by 
B — > X s 7 decay, which is already two order lower than the direct upper bound from CDF 
j^]. The result also implies that one may need data sample much larger than 10 fb _1 to hunt 
out t — )• cy signals at the LHC. 
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(a) (b) 

FIG. 5: (a) the Feynman rules of t^c interactions in the Lagrangian of Eq. [TJ (b) penguin diagram 
contribution to b — > S7 with top quark anomalous interactions. 



Appendix A: The calculation of C^{nw) 



Using the Feynman rules in Fig 5(a), the amplitude of penguin diagram in Fig 5(b) can 
be written as, 

iM = u s (p')[eT v (p,k)]u b (p)e u (k), (Al) 

r (p, k) = — r v cs v tb 



A cs 10 J (2tt) 4 [{pf - q) 2 - m 2 c + ie][(p - q) 2 - m 2 + ie][q 2 - m 2 w + ze] ' 

(A2) 

N = la L(j/f - </+ m q )a^(Kl R R + k^L)^- </+ m^Lg^k,, (A3) 
with R=(l + 7 5 )/2 and L = (1 - 7 5 )/2. By Dirac algebra 

la L^ u « R R + KLQifypL = 0, (A4) 
the terms with q 2 in N vanishes and N becomes 

N = m c Kl L W- q)a^ + (4 - D)a^(j/- ^}Lk^ 

+ m t Kl R [2a^(]/r - <f) + (4 - D)ty - q)a^]Rk^ (A5) 

Thus, there is no divergence in T u (p,k). After integrating out q in the T u (p,k) and using 
on-shell condition, T u (p, k) can be written in the following form, 

eV(p, k) = ie—^-y^ K%(mJ L (x)L + m b f R (x)R)] , (A6) 
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where 



K 



tcL 



A 



2m r 



(x c - l)(x t - 1) (x c - l) 2 (x c - x t ) 



^2m t 
A 



(x c - l)(x t - 1) (x c - l) 2 (x c - x t ) 



In x c + 



In x c + 



.77 



(x t - l) 2 {x c - X t ) 



X7 



(x t - l) 2 (x c - X t ) 



In x t 

(A7) 
In Xt 

(A8) 



Using the convention of Ref . 2l| , we have 



< )NP (M, 



1 V* s 
V£ * A 



+ 



_(rc c - l)(a; t - 1) (x c - l) 2 (x c - x t ) 



lnx,, — 



.17 



(x t - l) 2 (x c - x t 



In x 4 



(A9) 



Appendix B: Main formulas and inputs 



Following the notation in Ref. 20], the branching ratio of B — > X s 7 can be expressed as 



B[B -> X sl ] Ei>Eo = B cxp [B X c ez>] 
where P(E ) is the perturbative ratio 



v; 



r6 



6a P 



7T 



C 



[P(E ) + N(E )} 



(Bl) 



6a r 



\Vcb/V ub \ 2 T[B ^ X u eu] V cb n ' y un (B2) 
which includes the Wilson coefficients of Eq. [JJ N(Eq) denotes the non-perturbative correc 
tions. The semileptonic phase space factor 



C 



V 



ub 



T[B ^X c ev] 



(B3) 



Vd, T[B^ X u eu] 

can be obtained from a fit of the experimental spectrum of the B — > X c lv [23 [_. 

For calculating B{t — > c~f), we use the NLO formulas in Ref. [3] and 24j. Because 
t — > bW is the dominant top quark decay mode, the branching ratio of t — > cry is defined as 

T(t ->■ cry) 



^^ C7) -F(t->^)- 
The partial width T(t — >■ 07) at the NLO can be found in Ref. [if]], namely, 



r N Lo(£ ->■ C7) 



2a, 
9tt 



r (t ->• c 7 ) 



-31o,(^l-2,N 



(B4) 



(B5) 
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TABLE II: Experimental inputs for calculating the branching ratio of B — > X s j and t — > cy. 

Experimental Inputs 
a em = 1/137.036 [25J M z = 91.1876 ± 0.0021 GeV [25] 

a s {M z ) = 0.1184 ± 0.0007 [25] M w = 80.399 ± 0.023 GeV [25] 

G F = 1.16637 x 10- 5 GeV" 2 [25] mf = 4.67^;^ GeV [25J 

A = 0.8121°;°^ [26] m c (m c ) = (1.224 ± 0.017 ± 0.054) GeV [27] 

A = 0.22543 ± 0.00077 [26] m iiPoie = 172.0 ± 0.9 ± 1.3 GeV [25] 

p = 0.144 ± 0.025 [26] B cxp [B X c ev] = (10.64 ± 0.17 ± 0.06)% [28] 

fj = 0.342±g;gi| [26] C = 0.580 ± 0.016 [23] 

|V£ WK*| 2 = 0.9625 e ew = 0.0071 [20, 29] 

(Vus V ub)/(y t * s V tb ) = -0.007 + 0.018i N(Eq) = 0.0036 ± 0.0006 [20] 

V^/Vg = -24.023 - 0.432i E = 1.6 GeV 



where T (t — > cy) = amf (k^ cR / A) 2 is the LO partial decay width. 

The partial width of t — > bW has been calculated in Ref. j^J at the NLO, which reads 

2a „ ' 



r NL o(t ->• &W0 = r (t ->bw){i + 



3tt 



9-4/3 



ln/3^ + 2Li 2 (/3 2 



r / (i-^)(2^-i)(^-2) \ 2 

V ^(3 - W J ln(1 - ^ 

- 2Li 2 (l - 0*,) 



2/3^(3-2/3^) 



(B6) 



with T (t HP) = |^|^| 2 /3^(3 - 2/3^) and /3 W = (1 - m^/m 2 ) 1 / 2 . 

The experimental inputs are collected in Table. [TTJ in which the CKM factors are derived 
from the Wolfenstein parameters A, A, p and fj. 
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